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Using lattice QCD with dynamical Wilson fermions, we study 7 = 0 and J p = 0 + channel which is constructed 
by + d<7)|0), in order to search for the <r meson. Our preliminary result shows that the connected and 

disconnected diagrams contribute to the cr meson propagator in the same order. 


1. INTRODUCTION 

In QCD, the chiral symmetry is spontaneously 
broken (restored) at the confinement (deconfine¬ 
ment) phase where sigma meson plays an essen¬ 
tial role. We do not doubt this mechanism, and 
yet no one is sure whether the sigma meson ex¬ 
ists or not. It might have large mass and/or wide 
width, or might be simply 7r — 7t correlation effect. 
The light cr meson had disappeared from the ta¬ 
bles of Particle Data Group (PDG) for over 20 
years. However, the 1 = 0 and J PC = 0 ++ me¬ 
son, l ‘/o(400-1200) or cr”, appeared bellow 1 GeV 
mass region in PDG recently|jj. This is probably 
because of recent 7r-7r scattering phase shift re¬ 
analyses; Especially, Igi and Hikasa constructed a 
general model-independent framework, which re¬ 
spects the analyticity, unitarity and crossing sym¬ 
metry together with chiral symmetry low energy 
theorem, to describe the 7T7 t elastic scattering be¬ 
low 1 GeV mass region and investigated the ex¬ 
istence of cr meson §■ See Ref.|| for a good 
review on the situation and physical meaning of 
the sigma meson. 

Now it is very desirable to investigate whether 
cr meson appears as a pole based on lattice QCD. 
Using the quenched approximation, Alfold and 
Jaffe discussed the possibility of the light scalar 
mesons as q 2 q 2 states rather than qq ||. Mc- 
Neile and Michael computed the mixed iso-singlet 
scalar masses of qq and glueball states in two kind 


of situation, i.e, with and without the dynami¬ 
cal quark effects Q. The cr meson masses which 
are obtained with consideration of the dynamical 
quark effects are much lower than the quenched 
results. The goal of our project is to conclude 
whether the cr meson exists or not below 1 GeV 
in QCD. 

2. cr PROPAGATOR 

We construct 1 = 0 scalar channel by cr|0). The 
operator cr is given as 

a(x) = V L--■ « 

c—1 a =1 

where u and d are the w-quark and d-quark Dirac 
spinors, respectively. This operator has the same 
quantum number as the vacuum, 1 = 0 and 
J p = 0 + . The indices c and a denote color and 
Dirac spinor indices, respectively. The cr meson 
propagator is given by, 

G(y,x) = -(TrW~ 1 (x,y)W~ 1 (y,x)) 
+2(TrW- 1 (y,y)TrW~ 1 (x,x)) 
-2(TrW- 1 ( 2 /,y))(TrW- 1 (x,x)). (2) 

Here ”Tr” represents summation over color and 
Dirac spinor indices and W~ x is u, d quark prop¬ 
agator. The third term of Eq. ©, (<x(y)){<x(x)), 
corresponds to the subtraction of vacuum contri- 
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Figure 1. TrW~ 1 (x,x) evaluated by the noise 
method as a function of the number of Z2 noise 
sources. 


bution.Q The second and the third terms are the 
same order and the high precision numerical sim¬ 
ulations and careful analyses are required. 

3. NUMERICAL SIMULATIONS 

We calculate the a propagator by using Hy¬ 
brid Monte Carlo algorithm. We use Z 2 noise 
method to calculate the disconnected diagrams. 
The 1000 random Z 2 numbers are generated. The 
two-flavors Wilson fermion is simulated on the 
8 3 x 16 lattice. Based on ref. IQ, we set (3 = 4.8, 
k = 0.1846 (o = 0.197(2) fm, k c = 0.19286(14)) 
ji]]. After the thermalization trajectories, cr prop¬ 
agators are calculated on a configuration in every 
10 trajectories. 

3.1. NUMERICAL ACCURACY 

Since there are big numerical cancellation in 
(a(y)a(x)) — (a(y))(a(x)), we must be careful to 
controle the numerical accuracy. 

In Fig.[I], we show the values of a(x) for a 
typical configuration as a function of the num¬ 
ber of Z2 noise sources. Dotted lines represent 

1 (cr) is nothing but (ibip ). 


Figure 2. The connected and disconnected prop¬ 
agator parts of cr meson. 


the required accuracy to evaluate (a(y)a(x)) — 
(a(y))(a(x)). From the figure, we may conclude 
that 1000 noise sources are very safe and we do 
not suffer from any systematic error due to the 
Z2 noise. 

We have checked also the effect of CG solver 
tollerance, and T rW~ x {x, x) is very stable for e < 
HU 10 . 

3.2. PROPAGATORS 

As preliminary results, Fig| shows the con¬ 
tribution of the connected and disconnected di¬ 
agrams of the cr propagator in the time direction. 
The connected diagram shows clear exponential 
damping. On the other hand, the disconnected 
diagram has large error bars. However there ex¬ 
ists the contribution of disconnected diagram to 
cr meson propagator and its order is the same as 
connected diagram. Furthermore it can be possi¬ 
ble that the existence of the disconnected diagram 
makes the a meson mass lighter. This means that 
the estimation of the disconnected diagram is im¬ 
portant to determine the property of <r meson. In 
order to argue cr meson in detail, the improvement 
of the evaluation of the disconnected diagram and 
more statistics are indispensable. 

In Fig.|| we compare the cr meson propagetor 
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Figure 3. The comparison between the a propa¬ 
gator and p propagator. 


with p propagator. We can see that the a meson 
mass could be the same order of p meson mass, 
i.e., we obtain the suggestion of the existence of 
light a meson, though the error bar of a meson 
propagator is very large. 

4. CONCLUDING REMARKS 

We investigated the property of / = 0 and 
J PC = 0 ++ scalar meson (a meson) whose oper¬ 
ator is + dd). In the a meson propagator 

the contribution of disconnected diagram is the 
same order of connected diagram; Quenched ap¬ 
proximation is not reliable for the investigation 
of the c7 meson. The evaluation of the discon¬ 
nected diagram is done by using Zi noise method. 
A statistical error of a propagator which comes 
from the disconnected diagram mainly is large 
in the present stage. As preliminary results, we 
obtain the following properties of a propagator: 
(1) Both the connected and disconnected parts 
equally contribute to the cr propagator. (2) a 
meson could have mass of the same order of the 
p meson. 

It is necessary to generate much more gauge 
configurations and improve the statistical preci¬ 
sion of the estimation of a propagator. We also 


plan to improve the source of the sigma propaga¬ 
tors. 

Furthermore we must investigate the mixing 
state of the er meson and glueball if we obtain 
the result that er meson mass is greater than 1 
GeV regiong,?]. 
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